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Abstract 

The influence of hydrodynamic interactions on lane formation of oppositely charged driven col- 
loidal suspensions is investigated using Brownian dynamics computer simulations performed on the 
Rotne-Prager level of the mobility tensor. Two cases are considered, namely sedimentation and 
electrophoresis. In the latter case the Oseen contribution to the mobility tensor is screened due to 
the opposite motion of counterions. The simulation results are compared to that resulting from 
simple Brownian dynamics where hydrodynamic interactions are neglected. For sedimentation, we 
find that hydrodynamic interactions strongly disfavor laning. In the steady-state of lanes, a macro- 
scopic phase separation of lanes is observed. This is in marked contrast to the simple Brownian 
case where a finite size of lanes was obtained in the steady-state. For strong Coulomb interactions 
between the colloidal particles a lateral square lattice of oppositely driven lanes is stable similar 
to the simple Brownian dynamics. In an electric field, on the other hand, the behavior is found in 
qualitative and quantitative accordance with the case of neglected hydrodynamics. 

PACS numbers: 82.70.Dd, 61.20.Ja, 64.70.Dv, 05.70.Ln 
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INTRODUCTION 



The dynamics of colloidal particles dispersed in a fluid solvent is quite different from 
the ballistic motion of molecular systems which is described by Newton's law 
The viscous solvent both damps the motion of a colloidal particle and leads to kicks of the 
solvent molecules with the colloidal particle leading to Brownian motion if the time-scales 
of the molecular solvent is much faster than that of the diffusive motion of the colloidal 
particles. In concentrated suspensions, a dragged colloidal particle influences the motion 
of another particles via the solvent flow field. These so-called hydrodynamic interaction 
is typically long-ranged. While it can be neglected in colloidal suspensions of very small 
volume fractions, it induces significant corrections in the equilibrium and nonequilibrium 
dynamics of colloidal suspensions 0, [|5]. The equilibrium structures and phase boundaries, 
on the other hand, are unaffected by hydrodynamic interactions. 



In the past years a simple nonequilibrium phase transition 



B, q y, is 



9| has been dis- 



cussed in a binary mixture of colloidal particles which are driven by a constant external field 



id, llll, Q, [l3( . The drive is different on the two particle species and could arise from gravity 
and from an external electric field in the case of charged colloidal suspensions. Brownian 
dynamics computer simulations with neglected hydrodynamic interactions strongly support 
the scenario that - as a function of the driving strength - the mixtures undergoes a transi- 
tion from a mixed steady-state with anisotropic correlations towards a steady-state where 



macroscopic lanes are formed. 



mixtures 



0, fl Q, Q, 
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re transitions has been found for oppositely driven repulsive 



161 ] in two and three spatial dimensions and it seems to be 



a first-order nonequilibrium transition with a significant hysteresis in an order parameter 
which detects laning 10 L The general scenario occurs also in pedestrian dynamics 0, Q 
and in granular matter 2, . 

Recently, lane formation was observed in real-space experiments by Leunissen et al j^ . 



Equimolar mixtures o 
nary ionic crystals 22 , 



oppositely charged colloidal particles were prepared which form bi- 



23] . These crystals were exposed to a strong external electric field and 



the dynamics of lane formation was watched by confocal microscopy. Subsequently extensive 
3rownian dynamics simulations were carried out to map the nonequilibrium phase diagram 



241 ] . These simulations assumed a Yukawa pair interaction between the particles and ne- 



glected hydrodynamic interactions completely. A wealth of different steady-state structures 
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was detected. In particular different orderings were found in the plane perpendicular to 
the drive including a square, triangular or rhombic crystalline lattice of lanes, a network 
structure with a finite structural length and intermediate chain formation of lanes {24]. A 
rough estimate of the experimental parameters used in Ref. 2JJ reveals that the laned state 
observed experimentally indeed falls into the region where laning is expected to occur by 
Brownian dynamics computer simulations. 

In this paper we address the question how hydrodynamic interactions influence the sce- 
nario and the steady-state diagram of lane formation. Our motivation to do so is twofold: 
First, the experiments, of course, contain hydrodynamic interactions in their full glory, and 
therefore an inclusion of hydrodynamic interactions is needed for a quantitative comparison. 
Second, there is a principal need to understand in which direction hydrodynamics influence 
lane formation. In particular it is known 25] that the leading long-ranged term in the mo- 



bility pair tensor is screened if an electric field is applied since a charged colloidal particles 
is surrounded by counterions of opposite charge. That makes the hydrodynamic interaction 
significantly different from e.g. sedimentation induced by different buoyant masses of oppo- 
sitely charged colloids where the action of gravity on the microions can safely be neglected. 
In the latter case, the leading part in the mobility tensor at large interparticle separation is 
the unscreened Oseen tensor. It would be interesting to explore how far the steady-state is 
affected by hydrodynamic interactions in both cases of sedimentation and electrophoresis. 

We use Brownian dynamic computer simulations and include hydrodynamic interaction 
by using mobility tensors on the Rotne-Prager level. Both cases of sedimentation and elec- 
trophoresis are studied separately with an unscreened resp. screened version of the mobility 
tensor. The steady-state phase diagrams and the drift velocity are simulated. The simulation 
data are compared to that obtained by simple Brownian dynamics where hydrodynamic in- 
teractions are neglected. For sedimentation, we find that hydrodynamic interactions strongly 
disfavor laning. In the steady-state of lanes, a macroscopic phase separation of lanes is ob- 
served, i.e. the sickness of the lanes are of the system size. This is in marked contrast to the 
simple Brownian case where a finite size of lanes was obtained in the steady-state. For strong 
Coulomb interactions between the colloidal particles a lateral square lattice of oppositely 
driven lanes is stable similar to the simple Brownian dynamics. In an electric field, on the 
other hand, the behavior is found in qualitative and quantitative accordance with the case 
of neglected hydrodynamics. All lateral structures are reproduced and the topology of the 
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steady-state phase diagram is unchanged. 

The paper is organized as follows: in section II, we describe our model and simulation 
scheme for sedimentation and electrophoresis. Brownian dynamics simulation results are 
presented in section III. We conclude in section IV. 



THE MODEL 



We perform Brownian dynamic simulations to study an equimolar binary mixture of 
2N = 1024 oppositely charged colloidal particles of diameter a dissolved in a solvent fluid 
of shear viscosity 77 at temperature T and volume fraction = 2Nna 3 /6/ 3 exposed to an 
external driving field, where I is the dimension of a cubic simulation box having periodic 
boundary conditions. Henceforth, a serves as the unit of length and k^T, the thermal energy, 



as the energy unit of the system. To mimic the experiments by Leunissen et al [2ll . |26| the 
particles interact with an effective screened Coulomb potential (or Yukawa potential) plus 
a steric repulsion 14: 

Z 7 pi-^inj/o—i)) 
V ™ = V \l + 'Jw r iit , + « 

with Vo = 50A;bT the strength of the interaction potential and Zi = ±1 the sign of the charge 
of particle i. = |r< — Tj\ denotes the distance between particle i and j, where Tj is the 
coordinate vector. The inverse screening length k governs the range of the interaction and 
is determined by the salt concentration of the solution. The steric repulsion between the 
particles, that prevents the system from collapsing, is approximated by a repulsive (shifted 
and truncated) Lennard- Jones potential 



V h (n 



12 / n6 ^ 



if m < 2 1 / f V 

(2) 



else, 



with e = 4Vo/(1 + ko"/2) 2 . The constant external driving field that acts in opposite directions 
on the two different particle species reads as 



next 



Zife z , (3) 



where e z is the unit vector along the z direction and / = lbOk-^T/a is the strength of 
the external force. The external force is supposed to stem from either an electric field or 



a gravitational field accompanied with different buoyant masses of the oppositely charged 
particles. Though the external force Eq. ([3]) may in both cases be identical, if the charges 
and/or buoyant masses are chosen accordingly, the hydrodynamic interactions are not. 

The algorithm used to simulate the diffusive Brownian motion of the colloidal particles 
was proposed by Ermak and McCammon [27|]. Here, the translational displacements of the 
particles are deemed to occur in time steps of fixed length At and the update algorithm is 
given by j^ : 

N 



Ti(t + At) = r 4 (t) + At jr |Bsi|) . F ,(t) + V rj • Dji (t) X + Aif, 

j=l v B 



(4) 



where denotes the diffusion tensor field depending on the positions of the particles at 
time t. The random displacements Arf are chosen from a joint Gaussian distribution with 
mean and covariant matrix |28| 

(Aif > G = 0; (Aif Aif > G = 2D ij At, (5) 

where (. . .)g denotes the average over the Gaussian noise distribution. F«(t), % = 1, ...,N, 
comprises the nonhydro dynamic forces due to interparticle interactions, determined by the 
gradient of the interaction potentials in Eq. (jTj) and Eq. (|2J), and the external force Fff Eq. 
acting onto particle i. 

Hydrodynamic interactions are included in the simulation through the mobility tensor 
fJLij = Dij/kftT. In a first approach we neglect hydrodynamic interactions completely to 
asses its effect on the system. In that case the diffusion tensor is given by Stoke's law in 
diagonal form 

'Yfjiij = 5ijl, (6) 

with friction 7 = 37r^o"H, where <th is the hydrodynamic diameter. In the sedimentation 
and electrophoresis situation we approximate the mobility tensor by two-body interactions. 
In this approximation the divergence in Eq. (j4]) vanishes always [291]. When studying the 
sedimentation, the buoyant masses of the oppositely charged particles are supposed to be 
such that the same force acts on the two species but in opposite directions. The action of 
gravity on the microions can safely be neglected. Therefore, in sum no net force is acting 
on the solvent and overall it remains quiescent. Then, for a pair of spheres of hydrodynamic 
diameter <th the mobility tensor is approximated by the well-known Rotne-Prager expression 
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Sal + (1-6, 



'j > 



3cr H , 



fa) + 
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where 



A-(l + r®r); Q(r) = -^(1 - 3r <g> r) 



(7) 



(8) 



|r| 

with the unit vector r = r/|r|, ® a dyadic product, and 6ij Kronecker's symbol. On this 
level of approximation we incorporate all interactions up to 0((<7h/V) 3 ). Higher order con- 
tributions such as many-body, coupling between rotational and translational motions, and 
lubrication forces are neglected. The leading term in Eq. (jTJ) is given by O(r) which is of 
the order of l/|r| for large distances. 

However, when regarding the electrophoresis the mobility tensor has to be altered since 
forces induced by the surrounding counterions into the solvent are in sum equal to the 
force induced by a colloidal particle. Thus, the solvent flow stemming from the drag on 
the counterions cannot be neglected as is done in the sedimentation case. It results in an 
effective screening of the leading far-distance term of the hydrodynamic interactions between 
the colloidal particles as Long and Ajdai have shown [251]. Their mobility tensor /zg A reads 
as 



^1 + ^(1-5, 



1 a 



1 + — + — 



1 1 
— + 



2,. 2 I ;,) '" ; ' 



K 2 



(9) 



Here, the leading order term is Q(r^) which decays as l/|r 3 |. To account for the polymer 
coating on the colloidal particles that gives rise to the steric repulsion we choose <th = 0.9a 
throughout this paper. 

More sophisticated simulation techniques for spherical particles in an unbounded space 
including lubrication approximation for particles in close proximity and multipolar expan- 



sion methods are available 



31 



32 



33 



34 



351 ] . However, in the electrophoresis where the 



hydrodynamic interactions of the counterions become important explicit simulations of all 
colloidal particles and their counterions - 110 per colloidal particle in the experiments by 
Leunissen et al - are still beyond computational means. Therefore, we adopted the calcula- 
tions of Long and Ajdari to our simulations and compare it to the sedimentation problem 
on the same level of accuracy, i.e. the Rotne-Prager level. To our best knowledge this is the 
first Brownian dynamic simulations with the Long and Ajdari mobility term [9) 



Both mobility tensors [7] and [9] are long-ranged and thus require an Ewald-like summation 
in simulations analogous to Coulomb and dipole-dipole interactions. Details on the sum- 
mation and discussions about appropriate boundary conditions to the system can be found 



elsewhere 



36 



33, 



38 



391 ] . We applied the scheme suggested by Beenakker [39] and adapted 



it for /x^ A accordingly. The square root of the diffusion tensor, needed when calculating the 
random displacements in Eq. (jSJ), are obtained from a Cholesky decomposition: 



D = L L 



(10) 



where L is a lower triangular matrix and L T is its transpose. A suitable time scale for 
our system is tb = r )o 1 jh^T. The equations of motion including the external field are 
numerically solved using a finite time step At = 2 ■ 10~ 5 tb in all simulations. Statistics 
were gathered after an initial relaxation period of 20tb- The starting configuration of all 
simulations was a homogeneous mixture. 



RESULTS 



Order parameter and steady-state phase diagrams 

To assess the effect of hydrodynamic interactions on the lane behavior of oppositely 
charged colloidal particles we study a set of volume fractions and inverse screening lengths 
k* = kg and map out nonequilibrium steady-state phase diagrams for all three situations: 
hydrodynamic interactions neglected (A), electrophoresis (B), and sedimentation (C). 

A state of lane is thereby identified by a laning order parameter that is defined through 

where the angular brackets (...)t denote a time average. The local order parameter $j = 
[n\ — n ) 2 /(n\ + n ) 2 is assigned to every particle i, where the numbers n\ and n Q are 
the number of like charged particles and oppositely charged particles, respectively, whose 
projections of distance onto the plane perpendicular to the field are smaller than a suitable 
cut-off length scale z c . is equal to 1 if all particles within this distance criterion are of 
the same kind and zero if rt\ = n Q , i.e. a homogeneous mixture. We chose for convenience 
z c = |cx to detect all lanes starting from a single queue of particles. In what follows we will 
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FIG. 1: (Color online) Nonequilibrium steady-state phase diagram for a constant driving force 
of strength / = 150k^T/a with hydrodynamic interactions neglected accompanied by a typical 
simulation snapshot of the projection of the particle coordinates onto the plane perpendicular to 
the driving field for each different state. The lines between the phases are a guide for the eye. 

use a threshold: for $ > 1/2 we call the configuration a state of lanes while in the opposite 
case ($ < 1/2) we call it a state without lanes. 

We observe that lanes form different structures in the plane perpendicular to the driving 
direction for different values of k* and 0. We find lanes placed on a square or triangular 
lattice, a network-like structure (reminiscent of a bicontinous microemulsion or microphase- 
separated system), coexistence regimes of the same, and macroscopically separated lanes. 
The resulting nonequilibrium steady-state phase diagrams are shown in Figures [I], [2} and 
|3j They are accompanied with typical simulation snapshots of the projection of all particle 
coordinates onto the xy-pleme of the respective situation. 

What can be seen at first sight is that the qualitative behavior of situation (A) and (B) 
in Figures [1] and [2] is almost identical with only subtle differences, while on the other hand 
the phase behavior changes drastically for situation (C), Figure El In the latter the whole 
phase diagram is altered and the diversity of phases found is reduced compared to the first 
two cases. 
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FIG. 2: (Color online) Same as Fig. Q] but for electrophoresis with hydrodynamic interactions 
taken into account through jLt^ in Eq. ([9]). The phase diagram reveals only minor differences as 
compared to the the case of neglected hydrodynamic interactions in Fig. (TJ 




FIG. 3: (Color online) Same as Fig. [T]but for sedimentation with hydrodynamic interactions taken 
into account through /x?- p in Eq. ([7]). The phase diagram shows significant changes as compared 
to Fig. CD and Fig. H 



Comparison of simulation results for neglected hydrodynamic interactions and elec- 
trophoresis 

In this subsection we briefly describe the two phase diagrams in Figures [Hand [2] and their 
differences, beginning at low volume fractions and ending at high ones, and then dwell on the 
third diagram, Fig. [3J thereafter in subsection . A more ample discussion on how different 
phases are identified and what structural correlations they exhibit can be found in a previous 
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work of the authors on the same system with hydrodynamic interactions neglected but for 



a slightly different driving strength and larger systems [24|. Here, we find for situation (A) 
virtually the same results as in the previous work with only one difference, namely that we 
do not encounter a rhombic phase for = 0.4 and n* = 1, 2, 3. 

For very low volume fraction, <p < 0.01, in both systems the correlations between the 
particles are not sufficient to form lanes at all. Thus, the systems are in a phase of no- 
lanes. Only for very low salt concentration, i.e. small k*, where the electrostatic coupling 
between the colloidal particles is strong, we find a coexistence region between lanes and 
no-lanes. Here, the region with no-lanes consists of voids, where hardly any particle is 
found. The structure of the lane region, on the other hand, is different in the two situations. 
For situation (A) the corresponding snapshot in Fig. [T] reveals fixed lattice points while 
the snapshot in Fig. [2] situation (B) shows a network-like structure. For situation (A), 
an initial configuration with lanes placed on a square lattice separated from a completely 
depleted region is stable in simulations, as well. Thus, we assume that in situation (A) the 
lanes/no-lanes phase is a transient state toward a complete square lattice and no-lane phase 
separation. Hydrodynamic interactions destroy the coexistence phase for = 0.01. It only 
occurs in a denser system with <fi = 0.1 whereas in situation (A) this state shows already 
up at = 0.01. Additionally, the voids are more pronounced in the latter case compared to 
situation (B). Upon increasing n* > 2 for = 0.1 in both situations we find a network-like 
structure whose characteristic spacing is increasing with increasing k*. For situation (A) 
there is also a small coexistence region between network and square lattice at k* = 2. 

To obtain a quantitative measure of the characteristic spacing in the network structure 
we determine a structure factor perpendicular to the driving field of like-charged particles. 
The steady-state partial structure factor has been calculated by evaluating the expression 

S AA (k) = l+ph AA (k), (12) 

with p = 6<p/7ta 3 the number density and the wave vector k = |k|, where k = (2ir/l)(k x , k y ) 
and k x , k y are integers. h AA (k) is the Fourier transform of the total correlation function 
^AA(rj_) = #aa(?"±) - 1 with r = (rj_, z) and 

gAA(r ± ) = — ( Yl ^(r ± -\r l± -Tj ± \)5(zi- Zj)) , (13) 
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FIG. 4: Partial structure factor S±(k) of like-charged particles perpendicular to the driving field for 



k* = 10 and <j) = 0.1 for hydrodynamics neglected and electrophoresis. The pre-peak at k$ = k^a 



indicates an additional length scale of the structure in the network-like phase. The inset shows 



the position k$ of the pre-peak as a function of the inverse screening length k* for a fixed volume 



fraction of <j> = 0.1. 



where <5(x) denotes Dirac's delta-distribution. An example of the steady-state partial struc- 
ture factors for k* = 10 and <fi = 0.1 for both situation, (A) and (B), is shown in Fig. HJ One 
clearly observes a pronounced pre-peak at the wave number ko in both cases. A pre-peak 
in the structure factor is an indication of an additional mesoscopic length scale as typical 



for bicontinous networks, such as e.g. microemulsions [40L 1411] . In the inset we additionally 
present the position fc of the pre-peak as a function of the inverse screening length. It is 
evident from the picture that the characteristic spacing is indeed growing with increasing 
K*. We find hardly any difference between situation (A) and (B). 

For = 0.2 an additional phase for small inverse screening length shows up in both 
phase diagrams. Oppositely driven lanes are placed on a square lattice with an alternating 
charge pattern. The formation of this lattice structure can be qualitatively understood 
from an effective interaction between oppositely charged driven lanes which has a short- 
ranged repulsive and a long-ranged attractive interaction. The former stems from the friction 
between oppositely driven particles while the later results form the Coulomb interaction. 
The square lattice then reduces the electrostatic energy of the system because each particle 
has only oppositely charged neighbors. For increasing salt concentrations we encounter a 
coexistence region between the square lattice and the network-like phase and finally end up 
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in a pure network-like phase. The phase diagram is in both situations very similar, only the 
borders of the transitions are slightly shifted. In the electrophoresis case the network-like 
structure is preferred to the square lattice. 

For a higher volume fractions of <fi = 0.3 a coexistence regime between a triangular lattice 
and a network-like structure is found. The lattice-points in the triangular phase are rather 
randomly decorated with different charges. Here, the short range repulsion plays the dom- 
inant role compared to the electrostatic interaction. It enforces a triangular lattice due to 
packing effects although electrostatically it is strongly disfavored because like-charged parti- 
cles necessarily occupy lattice points next to each other. Again, hydrodynamic interactions 
slightly shift the phase boundaries towards the network-like structure. 

For the highest volume fraction studied, <fi = 0.4, both phase diagrams show exactly the 
same behavior. Here, the short range repulsions dictates the phase behavior for nearly all 
salt concentration but for k* = 1 and enforces lanes to be placed on a triangular lattice. 
Only for k* = 1, where electrostatic interactions are prominent, a square lattice is preferred. 
In principle a square lattice is possible up to the packing of a simple cubic lattice of <fi = 0.52. 

In summary we observe very similar behavior in both situations. The observed differences 
can be qualitatively explained by the fact that hydrodynamic interactions disfavor lanes 
driven oppositely past each other. 

Sedimentation 

Regarding sedimentation, Fig. [3l the whole phase diagram exhibits only three different 
phases. For volume fractions < 0.1 we do not find lane formation for all inverse screen- 
ing length studied. For increasing volume fractions and strong electrostatic interactions, 
K* < 2, first the square lattice at <p ~ 0.3, that is also present in the previous two situa- 
tions, is recovered and then the system reenters a region with no-lanes for = 0.4. This 
behavior nicely illustrates the competition between hydrodynamic interactions disfavoring 
lanes driven oppositely past each other and the electrostatic interactions favoring a square 
lattice. Only for the small regime around <fi ~ 0.3 the electrostatics succeeds the hydro- 
dynamic interactions and enforces a square lattice. For all other volume fractions laning 
is destroyed. However, for stronger salt concentrations, where the Coulombic coupling is 
reduced, we discover a situation that is not present in the previous situations (A) and (B), 



12 



140 - 




a □ electrophoresis 
o- O sedimentation 



80 - 



■ . 



2 



4 



5 



6 



7 



8 



9 



10 



K 



FIG. 5: Average dimensionless drift velocity v* = vt-q/ct in drive direction as a function of the 
inverse screening length k* at (ft = 0.3 for Brownian dynamic simulations with hydrodynamic 
interactions neglected, taken into account through /J^j, and nff ■ 

namely a region in which only two big completely separated lanes. We call this state phase 
separated. In that case the long ranged hydrodynamic interactions prescribe the structure 
and the short ranged Yukawa interaction plays its role only at the rough interface of the two 
phases. From our simulations we conclude the lanes are separated by half of the box length. 

Drift velocity 

Now, we study the influence of hydrodynamic interaction on the drift velocity along the 
field direction that is defined as follows 



This entity measures the mean-square displacement of each particle in the nonequilibrium 
steady-state. A study on the effect of hydrodynamics on the drift velocity of like charged 
colloidal particles was carried out by Watzlawek and Nagele [42]. We study two cases, first 
we fix the volume fraction at (ft = 0.3 and vary the inverse screening length and afterward 
vice versa for k* = 1. 

In Fig. Owe display v* = vr^/a for a fixed volume fraction (ft = 0.3 as a function of the 
inverse screening length k* for all three situations. For all cases the drift velocity increases 
with decreasing Coulomb coupling because oppositely charged colloids attract each other 
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FIG. 6: Average dimensionless drift velocity v* = vtq/ <j in drive direction as a function of the 
volume fractions <f) at k* = 1 for Brownian dynamic simulations with hydrodynamic interactions 
neglected, taken into account through M^ A , and nff- 

while driven in opposite direction and lanes mutually retard each other. For very strongly 
screened particles where this friction is less important all three curves reveal approximately 
the same value of v ~ 130tb/o~. Accordingly, this value is close to the drift velocity of 
v = 150tb/V for a system of infinite dilution subjected to the same driving force. While 
in (A) and (B) v grows gradually, in the sedimentation curve we encounter a jump in the 
drift velocity between k* = 3 and k* = 4. This coincides with the transition from the 
no-lane regime to the phase separated regime, see the phase diagram Fig. [31 On the other 
hand for n* < 2, where we find a square lattice, the drift velocity is similar to k* = 3. 
From that we conclude that the phase separated state of lanes supports particle transport 
while lanes placed on a square lattice enforced by strong Coulombic interactions slows down 
particle transportation. A further interesting feature is that curves for (A) and (B) intersect 
between k* = 1 and K* = 2 and that the screened hydrodynamic interaction enhance 
the drift velocity for larger inverse screening length. The same is true for the unscreened 
hydrodynamic interactions in the sedimentation for k* > 6. When studying the drift velocity 
for a fixed inverse screening length but for varying volume fractions in Fig. [6] we find again 
an intersection point of the curves for situations (A) and (B). Here, the drift velocity in 
the electrophoresis reaches an approximately constant value of v* ~ 112 for <p = 0.1 — 
0.4 whereas it gradually decreases when hydrodynamic interactions are neglected. For the 
sedimentation the drift velocity decreases monotonically. In contrast to the case of varying 



14 



salt concentration we do not encounter a jump in the sedimentation drift velocity when 
entering the square lattice at <fi = 0.3 and reentering the no-lane regime at <fi = 0.4. 



CONCLUSIONS 



In conclusion the influence of hydrodynamic interactions on lane formation of opposite 
charged colloids driven by an electric field or by gravity was investigated by Brownian 
dynamics computer simulations. Hydrodynamic interactions were included on the Rotne- 
Prager level. For an electric field, the leading Oseen term is screened due to the presence 
of counterions. The latter fact has lead to very similar steady state phase diagrams for an 
electric field as a driving source than that in the simple case of neglected hydrodynamic 
interactions. Various steady state were obtained as a function of the colloidal density and 
the range of the interaction. They can qualitatively be understood in terms of a competition 
of the mutual Coulomb attraction and friction of sliding lanes. At high densities the lateral 
structure is crystalline, the crystal is either triangular as dictated by packing at high densities 
and high screening or square-like at low-screening which minimizes the Coulomb attractive 
energy. On the other hand, in sedimentation where the two colloidal species have the same 
buoyant mass up to a relative sign, friction of sliding lanes is strongly enhanced leading to 
macroscopic separation of lanes. 

The steady-state phase diagram can in principle be verified in real-space experiments 
of charged suspensions which are driven in an electric field or sedimenting 21j. It would 
be interesting to construct a microscopic theory for the lane transitions which includes the 
lateral crystalline structure. The instability analysis within a dynamical density functional 
theory as applied to the case of equal charges in two spatial dimensions 12|, |43j should in 
principle be generalizable to the case of oppositely charged particles. 

Finally more sophisticated simulations schemes are needed in order to go beyond the 
Rotne-Prager level of approximation used in this paper. ^Among the various promising 
approaches are the stochastic rotation dynamics code 
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including hydrodynamics 



43, 



fluid particle dynamics methods 51 



48 



49 



52, 



451 ]. a lattice Boltzmann theory 



50l| and counterion flow or the recently developed 
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